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ABSTRACT 


Based on the theory of continued fractions, a technique 
is developed for the reduction of high order multivariable 
systems. The mathematical basis for which these techniques 
work is elucidated, and its superiority of the mixed form 
over any other form of continued fractions is established. 
The general solution to linear regulator problem is developed 


and the properties which this solution exhibit are elucidated. 
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Teen ODE Cl LON 


The capability of obtaining simplified mathematical 
models for use in the analysis of high order dynamic systems 
has traditionally relied on the experience and ingenuity of 
the analyst. Usually, these efforts have been achieved 
using both frequency and time domain techniques. 

In dealing with the problem stated above, and through 
the use of theory of continued fraction, this research has 
found a series of properties applicable to reduction of 
multivariable systems, lower order observers and derivations 
of lower order systems for the linear regulator problem. 

The relationships developed, show applications in areas 
where the high order systems are impractical or undesirable 
memmise due to their complexity or difficulty of implementation. 

In Chapter II reduced order models of multivariable 
systems for the first, second and mixed Cauer forms are 
developed. Techniques for approximating a high order linear 
time-invariant system with various inputs and various outputs 
by a reduced order model, have been suggested by Chen [1], 
Meier, L and Luenberger [2], L. S. Shieh and Y. J Wei [3], 
Meee. Calte and M. Healey lt], L. S. Shieh, J. M. Navarro 
EUER Yates [51,; D.—A. Wilson [61, L. S. Shieh and F. F. 


Esudjano [7]. 


cO 





Most of these methods for reducing high order linear 
systems are based on the following principles: 

1. The low performance terms can be discarded and the 

high performance terms should be retained. 

fee Limear transformation to obtain matrix diagonaliza- 
tion where certain diagonal elements can be dis- 
carded. 

3. The sum of squares of the errors between the 
responses of the real system and those of the 
approximate model at the sampling instant is 
minimized in order to obtain the parameters of the 
approximate model. 

In a recent paper, Chen [1] proposed a reduction of 
multivariable control systems by means of matrix second 
Cauer form of continued fractions. Through the method, a 
simplified model is obtained by keeping the first several 
significant matrix quotients and discarding the others. 
However, the technique (due to the nature of the Cauer 
second form), provides satisfactory results in the steady 
state region only. Furthermore, M. R. Calfe and M. Healey 
[4], have shown that the method does not guarantee the 
reduced model to be stable. 

In Chapter III derivation of lower order system for the 
linear regulator problem via Cauer form is obtained and 
also a near optimal solution for the original system can be 


found through a reduced system. 


m 





IA Outen us SYSTEMS REDUCTION VIA THE CAUER FORM 


Pee REE MATRIX CAUER FORMS 

L.S. Shieh and F.F. Gaudiano [7] have shown that in terms 
of multivariable systems the quotients in the three Cauer 
forms are replaced by matrix quotients and the division in 
the continued fraction is replaced by matrix inversion. The 


Biest matrix Cauer form is 
E pes E E P 
T(s) = (His + [HY + (His + [HL + [... III T IT (n 


mmremsecond matrix Cauer form is 


_ 1 rer ee 
Res) = CH, + [H,Z + CH, + s Eu]: rer 22) 


wehe mixed matrix Cauer form is 


i 
+ S f = 
LX, ee Kus 


= i ! 
ECOL EK. + Kis + [KZ +K > 


it 2 


+ KL + Bea Eee (3) 


where HT» Hs > Ky > and e are constant m x m matrix quotients 
obtained respectively from the matrix Routh's array and the 
generalized Routh's algorithm shown in equations (Ha) and 


(4b). 





Ha = A 
H, = A 
Ha = A 
0] B 
fees A 
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ae? Fo in eal 
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a >=: i SER pi 
ee, E 

31441 < >: 7 TE 
Aud Aug i E 

E a qe 
i rk ES HFA, ann 
ns 
BE u > ee 
en 
(4a) 
ES Ts en ee 

11821 < Du " SEMEL 
TOME A52 m 

21 = < > PL E 
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iA < b. y= Ag nA ne 
o er, 
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RN ES 
E -1 X E -1 
EE in? ENT < eal en 
en n 
K =A, 14+ DK =A „Æl 
ngon" pt. no 2m us 
on 
(4b) 


where the elements of the first and second rows of equations 
(4a) and (4b) are the matrix coefficients of the system 


given in equation (5). 


er A 


l |x 
2 23 22 gu 


T(s) = [A,,n Sr 


(5) 
n n-1 2 -1 
[A, ,n*18 + A, ons +...A,,5° + AS + A] 


The elements of the third, fourth and subsequent rows in (4a) 
and (4b) are evaluated respectively for the three Cauer forms 


"une formulation shown in (6a) through (Se). 


Hb = Ap,n+2-pAp +1,1 p 
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pa 
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ÉS 





Hp = Ap,lAp +1,1 Dan 


Fee An 2, 1810 - HA, -2A)-1, atl 3=3,45...2n (6b) 
det Ap+l,1 £z 0 N ERN 

I ES Ap. EET DD ELE E 
Kp = Ap, n+2 - p Apt, ntl-p Veen 


EE te Ke 1 ae 22 749-1. 1 i=1l,2,3, 


(6c) 


It is important to note that since the Cauer first and second 
forms are special cases of the Cauer mixed form, their 
formulation in (6a) and (6b) can be derived directly from 
(6c), by letting all kp's or all Kb's equal to zero respec- 


tively. 


ESTATE SPACE FORMULATION FOR THREE CAUER FORMS 
The Cauer Mixed Form - Consider a typical feedback system 
with a minor feed forward loop as shown in Figure 1. The 


overall transfer function is given by: 


mis) = is = LG) ae Libs (G 4 F)H] +. CF: 


Equation (7) can be rewritten as a mixed matrix of continued 
fractions 


meee nee, (8) 
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where 


= + 
H Ka KS 
D e (9) 
2S 2 
ne erdeneıey Matriz 


ICs) 





Figure 1. Feedback and Feedforward controls. 


If the subsystem G is expanded again, equation (3) is obtained. 
This equation can be represented by the block diagram shown 
in Figure 2. Where a 2 inputs-2 outputs nth order system 

is shown. Again, it is important to note that if all Ki's 
go to zero in Figure 2 the block diagram representation of 
Cauer matrix first form nth/2 order system as shown in 
eeure 3 will automatically be obtained. Ina similar 
fashion if the k'i's go to zero, the block diagram represen- 
tation of Cauer matrix second form nth/2 order system as 


Shown in Figure 4 is obtained. 


de» 








Figure 2. General Matrix representation of a nth order system by 
Cauer mixed form with two inputs two outputs. (n even) 


T6 








ystem by Cauer 
. (n even) 


Matrix representation of a nth/2 order s 
first form with two inputs two outputs 


Figure 3. 
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Figure 4. Matrix representation of a nth/2 order system by Cauer 
second form with two inputs two outputs. (n even) 
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Going back to Figure 2 and allowing ei to be the state 
variable vector (same order as the matrix K or K'), the 


equations in the time domain can be written as follows: 


E. DTI T E 
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Y = Kye, *Khe, tk e5tk, 2,7 Tr *K e 2 M A S2" (14) 
where U and Y are the input and output variables of the 
system respectively and I is the Identity Matrix. From 
Bouacıons (10), (11), (12), (13), and (14) the following 
corresponding matrix formulation can be obtained: 
BEA JE = FA-]E * EBJU, (15) 
and 
Be (01 es Dos. (16) 
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gin E ES 7 
an X(n/2)m 
(26) 


2 





| 
— 


mei *5m*l £(n/2-1)m+1 
= : ^m*2 *omt2 en? 
j sa a DS : 
^m “2m *3m |*(n/2)m 
(27) 
Ax ntl Xomtl *(n/2-1)m*1 
d^ NE» — | *2m*2 * (n/2-1)m*2 
B ee Pj ny . 
^m *^m * am *(n/2)m 
(28) 


where m is the number of inputs and outputs. It is important 
to emphasize the properties that the Cauer matrix mixed form 
exhibits. If again the Ki's go to zero in the state formula- 


mom described by equation (15) and (156), they will reduce to: 


EIE = -LA, , dE LOO 29) 


TL 


and 
E. (30) 
and letting K'i's - H'i's equation (29) and (30) define the 


state space formulation of the Cauer matrix first form for 


an nth/2 order system, wnere n is even, Equations (29) and 
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(30) are second order differential equations which can be 


simplified to first order differential equations by 


assigning a new state variable [2] 


(29) and (30) can be rewritten as: 


LA, y 12 : [A,, 12 + [BJU, 


and 


T 
fe Zu 
LC, d 


mS. hi. 
where LA, 1 1 Lom [A] and K'i's 


Kk uc 0 
then: 
1 t tI t t 1 
HoH, As Ay HSH, 
t t 1 t t 1 
HH, (CH, +H’, Hl, CH) +H', Hl, 
= 1 t t 1 t 1 1 
LA, 115 HH, CH! +H}, ) Hi, (H3 +H} +H} )H, 
t 1 1 t 1 1 1 
H-H, (HS HA 2H], CH) +H +H HE. . 
Ze 
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ns equations 


(31) 
(32) 
ptus 
1 
.HYHL 
1 t 1 
. CH +H! HE 
T! 1 TT! 
. (HI +HS+ALIH, 
t t p t 
CH) +H fH HL), 
(33) 
(34) 





where C42 877 € +: €, /9 “re given in equation (27). 
[A,, 1 = Nm En SEN (3/55) 
Jans) 


where I is the Identity matrix. 





Lr 
11702 
(36) 
where 94,5 95» €4 ... 94,5 are given in equation (28) 
De il =F eae eC te mas das US E A O HNIC 
21L 2 2E 2 4 ` ie) 
Ki's>0 (37) 


Mes tormulation given by equations (30) through (37) could 
have been obtained directly by inspection from Figure 3. 


BSmilarly, 11 all K':'s go to zero in equations (15) and 


(16) then 


EE -[A,„JE A (33) 
and 


MOMS E. (39) 


Zo 


Al » 
x 





where 


Ew pos limLA, J aadi sss a's 


32 
K'i's^0 
A, A, A, Ay H4 He co HH 
H-H, (H; *tH5)H, (H4 *H4)He OOHORFHAUH, 
E Is HH (H. *tH4)H, (H,*H 


+Ho JH... CH) +H 


2 3 


gtH, HL 


17 D 
AAN (A, +A, (H+H +H,.9H,...(H+H,+H ae )H 


9 3 5 aan 
L m 
(40) 
and 
Bee and eis 
1% 1 
K'3 s 
EC E = [E H, H vs (HI) 
12 pau p come 


which define the state space formulation for the Cauer matrix 
second form. Again, this preceding matrix formulation could 
have been derived directly by inspection from the block dia- 
gram shown in Figure 4. Equations (40) and (41) have the 
same form of the state space formulation given by Chen [1]. 
At this point, the following observations may be made re- 
garding the state space formulation just developed: 

1. The elements in the state matrix [Ai's] are simple 


matrix combinations of the matrix quotients 


2 











obtained from the continued fraction expansion or 
ROUN S algorithm, 

2. The submatrices appearing below the main diagonal 
have the same value as the submatrices at the 
diagonal (with exception of [A,]). 

3. The submatrices which appear above the main diagonal 
can be expressed in terms of matrix quotients in a 
very regular way. 

bh. State space formulations for the first and second 
Cauer form can be obtained directly from the Cauer 


Miedo nero ue direet Substitution. 


P PROXIMATION BY THREE CAUER FORMS 

The reduction of the order of a transfer function or 
decreasing the dimension of a state matrix is highly 
desirable or sometimes necessary in the analysis and design 
Seeecontrol systems. 

In terms of continued fractions, the simplification 
problem is carried out by expanding a given transfer function 
into one of the Three Matrix Cauer Forms of continued frac- 
tion and ignoring some matrix quotients. If the given 
system is outlined in state variable form, the simplifica- 
tion method is realized by partitioning the matrix and 
discarding some parts. Several examples are included 
for demonstrating the power of the method. Also a thoroush 
comparison among the three Cauer forms is presented and their 


advantages and disadvantages are discussed. 
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Feedback Gain and Feedforward Gain - Consider the system 
shown in Figure 5. The closed loop of the overall transfer 


Emmesıon is known as follows: 


-1,-1 





E 0 LES + Gl (42) 
Figure 5. Block diagram for a typical feedback system 
with a minor feedforward loop with two in- 
pues mama ewe COLT DUES: 
where G = 0. Equation (42) can be considered as the simplest 


continued matrix fraction expansion. The physical meaning 
implied in the formula is significant. It is easily seen 


that when E, is high the overall gain can be approximated by 


E 
l 


System. 


usnm other words, H4 dominates the behavior of the 


2g 





This fact establishes the feedback loop as the most im- 
portant link for influencing the behavior of the systen, 
leaving the feedforward loop as the second most important 
link. Furthermore if the subsystem G, instead of being zero 
Memstill a high order transfer function, it is possible to 
continue the expansion one after another. This corresponds 
to a combination of many feedback and feedforward blocks as 
shown in Figure 6. It should be noted that the most dominant 
term is A, and the second influence term is Pye When the 
matrix quotients in the continued fraction are lower and 
lower in positions, they are less and less important as far 
as the influence to the performance of the system is con- 
cerned. This observation is the general basis for the sim- 
plification technique developed for multivariable systems. 

Considering a simple case, a second order transfer 
function such as: 


s + A. 


= 2 
T(s) = [A,,S + A, JA, 35% + Aj95 + Aj, 


(43) 


can be expanded into three different matrix Cauer forms of 


Bentinued fraction as follows: 


F](8) = [H^ S + LH, - Dore * iie s (44) 
E 1 ee 

JP = LH. + Tal ~- CH, + E I (45) 
_ 1 Eu 

Be = LK, + K* S * [K,= + K J ] (46) 
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Block diagram representation of a continued fraction expansion. 


Figure 6. 


al 





where e Ca) ab ees haan F,(8) are the first, second and 


Hr EZ 


mixed Cauer continued matrix fraction expansions of equation 
n3). 
From previous observation, it is known that the most 


dominant term in equations (+4), (45), and (46) is His, H 


and Ky T Kis respectively. It 1s desirable to find a 


meaningful interpretation for the dominant term of each one 


1 


of the expansions just performed. This task is accomplished 
by applying the initial and final value theorems to equations 
un... (45), and (46) as follows. 
ie cavier First Form 
Performing the inverse procedure on F7]18) it can be 


Kpıtten as 


Zr at est pt 
FP (8)= 57735" 1H) HH) S+H)+H), J0H7 H HH) S * (HS HP FHA Hy *HAHI, 





)5*1] 1 
(47) 


Applying the final value theorem to equation (47) and 


allowing 


1 : ; 
U(s)= = where m is the number of inputs is 


Pound that 





m 








j Geo cuts ymo 
s>0 | 
dl 
pm 
(48) 


similarly, applying the initial value theorem and allowing 


| 
l 
J E 
A 
3 >00 
1 
m 


(49) 


The meaning implied in equations (48) and (49) is very 

significant. The initial conditions dominate the behavior 

of the system. In other words, the Cauer first form 

influences very heavily the transient part of the response. 
AEC auer second. Form 


Performing the inverse procedure of F2 69), this can 


Be written as: 
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mas). 2 ee Seal 
F968) "acy LH, tH, StH HH, IETS *(H4H,*H,H, *H4H, )StHHAH4H, ] 
C50) 


Applying the final value theorem to equation (50) and allowing 


pu 


C 

om 

0 

Du? 

i 
(Q [i 


it oe ound tchat y tss LA, 


s>0 i 
E 
"m = 


(51) 
Similarly, applying the initial value theorem and allowing 

5 

= Ei 

=? 22 

D= |’ y (t>0)=[H, + H ] 

: ^ $-o 4 

i > 
(92) 


Pi e-weesults obtained in equations (51) and (52) imply that 
the final or steady state value dominates the behavior of 
the system. In other words, the Cauer second form influences 
very heavily the steady state part of the system response. 

3. Cauer Mixed Form 


Performing the inverse procedure on Fag CS) this can 


be written as: 


MS TE i x Sl 
Fea CS? 7gts y LESS* KP JE KS ce eu ck K5)S*K. K,] (53) 





[EP 
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Applying the final and initial value theorems to equation 
(53) in the same fashion as for the previous forms, the 


following results are obtained: 


BE: lle ie 
y(t>») = [K] K 


s>0 : 
I 
m 


(54) 
and p 
E 
_ -1 
t0) = [K] 1, 
S>% 
a 
m 
(55) 


These results show the steady state value and the initial 
conditions to dominate on equal levels of significance the 
behavior of the system. Thus the Cauer mixed form influences 
the system response in the transient part as well as in the 
steady state part. This fact makes the Cauer mixed form a 
better and more accurate device to be used in the simplifi- 
cation and reduction techniques of transfer functions of 


multivariable systens. 
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Pree Oe UNCATION OF THE CAUER MIXED FORM THE BEST APPROXIMATION 

BIS NDEEDNIDS MATHEMATICAL SENSE FOR MULTIVARIABLE 

Soli REDUCTION 

im a recent paper, M.J. Goldman and C.T. Leondes [8] 
developed the mathematical basis for the sinpi cation 
technique involving the truncation of the mixed Cauer form 
for the single input single output case. In this section, 
an extension of their work to the multivariable case is 
presented and its superiority over any other form of con- 
tinued fractions is established. 


Brest ransfer function matrix TCs) for multivariable 


system can be expressed as follows: 


A n-1l n-2 2 
T(s) = LA) pS t Ag p-15 t000*A9348 + ALS + And X 


n-1 2 -1 
LA, n+1° + Ay n5 +... Ay 35 + A, 75 + A, 1) 


(56) 


where Ai,j are constant, m by m matrices and Aij = aj[I], 
]- 1,2,...ntl, where each aj is a coefficient of the 
common-denominator polynomial or 


+ 
ms) = 2 ajs)7? and [1] is the identity matrix. 


The nth convergent of a mixed matrix Cauer form can be 


represented by the following two configurations: 
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Pee =k +K'S + S[K.+KiS+S[...SLK «X 17 5] 1175174 (57) 
n 1 2 mon 


n ii 2 
E oun heu. del -14-14-1,-1 
en alk, Stk, tSLK SK tel... glk ie y Kl 17171 (58) 
where 


A - KA UI + 5A 
nn 


n n n-1 -1 neo) T 1l 
(59) 
= + i + = = + K! 
B K Ba SK B1 SB i2» B 0 BS Ky cao 
ee ont tat EE Ye 1 ot 
IEEE 505 uu s 
(60) 
mee 1 + KI Rt HS to - ea | 
en nr ae 7 
Recurrence relations (59) and (60) have been derived from 
standard results in the theory of continued fractions 
(Rice 1964) and Ks and = from the generalized matrix 
Routh Algorithm (Shieh and Gaudiano 1974). 
Since 
Dow cod y (61) 
nn Den 
then 
A E SON 
n n 
B, Sue ms ccs (62) 
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PROOF: 


Substituting equation (62) into (59) 


DN C A A A 


a 1 n-2 


and Co) 


are E . + Skis "+ 
n Mans n 


1 
A zm 


2 


Dividing equation (63) by S™, recurrence relation (60) is 
obtained. (It is important to notice the generality of 
Euusesons (59) and (650). By making all the K's = 0, the 
recurrence relationships for the nth convergents for the 
Cauer first form is obtained, similarly by making all the 
K'is = 0, the recurrence relationships for the nth covergents 


une Cauep second form are obtained.) 


From (59) 

Rn p e 

i y LA, ~ In ET SAn-21A2_1 
and (64) 

-1 
= pe 1 ES 

E ien SK B i SEE 
Solving and simplifying terms 

Ae a i -S(A _1Bn-2 E A (65) 
BO by induction 

DN NENNEN (29)0s 27. (66) 


n n-1 n n-1l 


where I is the Identity Matrix. 
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Hence, the difference between two consecutive convergents is 


given by 


-1 = 1 E n-1l -1 
A B. - A B = (-S) [B B. 14 : (67) 


Similarly, from (60), and following the same procedure 


l 


! nt! mont xt 2 oL(Oonx t CSI t 
Anl BA (a 1-2 Bel n-2?: (68) 
Eom induction 
aB -BA = lsn, (69) 
n n-i n n-1l S S 


and the difference between two consecutive convergents is 


given by 
=>] rl ICH, er me) -1 
SD d mec S ( m [B Bh- (70) 


> 


By looking at the recurrence relations (59) and (60), is 
observed that when S=0 or S=w, respectively, all B and all 
P ape non-zero provided that all A and N are non-zero. 
This is equivalent to apply the final value theorem to 
equation (67) for a unit step inputs; and applying the 


initial value theorem to equation (70) for an impulse 


response. 
Hence, 
d) tA mE NER Seen ee essa) sD Bae le 
eta n-l'n-1:$S-0 S MN 
s>0 
= 0 A o AZ 


Con 


99 





A 





Similarly, 


a) t rl t rl - UR m ik t t zn 
HA n DS MA 
i S 90 

= 0 OO Es , n-2 


(52220 


Equations (71) and (72) show the results that were expected, 


namely: 


The (n-l)th approximant and the nth approximant 


goes to zero in the sense of minimizing the follows semi- 


norms of the difference 


n- a 
ra EE O 0 
n-l n-l 320 Eh J n-l s=0 
n-2 qs | 
! = — F = 
Nr Gg l E | ded e D 
: a S 35 
where 
a di -1 
en een 
and (74) 
t - nm B — At Te 
Fh? z AB, = Ed 
Note: A semi-norm is a norm which does not satisfy the norm 
moron ECs) || = O implies F(s) = 0". 


Since the derivatives in (73) correspond 
cients in the Taylor series expansion of the 
F Coana MS about the Doints S0 and 

n-l mesh 


can be deduced that the output difference to 
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to the coeffi- 
funetiens 
ir 


S=w. Then, 


a nth approximant 








has the same Taylor series as the original function for 


n-1 


terms up to S ands. (=) ER 


= , effectively neglecting the part 


ene transfer function which differentiates the input > n 
times in the steady state as well as in the transient. This 
gives the sense in which the approximation by the Cauer mixed 


Torm works. 


BEE STEHPLIFYING A MATRIX TRANSFER FUNCTION 


If the following nth order system is given, 


E n-1 2 2 
T(s) = [Ay pS t ...An,S + A,,S + A,,S + A7q 1 X 
n 9 2 -1 
[A] n+1° + ...A,,8S° + AS. +A.,S + Aj] ”. 
(75) 
a simplified model of the system is desired. By performing 


the Generalized Routh's Matrix Array, the matrix quotients 
of the Cauer mixed form in equation (3) can be obtained. 

If an m order is desired only the first m pairs of K's 
and K'is should be kept in equation (3) and the remaining 
should be omitted. After the inverse procedure has been 
performed on the truncated continuted fraction, the simpli- 
fied model is obtained. 

For example the general transfer function obtained by the 
matrix continued fraction expansion for the Cauer mixed form 
1s given in equation (76). 


= 1 lie ' Le e 
AS E EK K SALK SK E 1] 3 3 1 > 


ul 








If a second order simplified model is required, only the 
first two pairs of matrix quotients will be kept, that is 
and the rest should be discarded. 


STE eae and K! 


MM" 2 
After the inverse procedure has been performed on the 
truncated matrix continued fraction expansion, the transfer 


function given in equation (77) is obtained. Equation (77) 


is the simplified model in Mixed Matrix Cauer Form. 


2 -1 
E t !oy' ! ! 
T(s) = [K,S+K,] x EEE t(I+K K +K K,)S+KK,] (77) 


The truncation procedure outlined above, applies in similar 
manner for the Cauer's first and Cauer's second forms. 

This methodology is particularly advantageous when state 
Space terminology is used. For the Cauer mixed form, the 


state space formulation is written as follows: 


wa KK AS AS 





ie 14 1^6 In 

1 t T 1 1 f 1 1 1 1 
ESOS! FRIKI TINO KURT 

ct ut 1 f ! ! f ! f f ! pct x £ 
OK (KS KA KE... (KD +KE+KL DK E, 
f T wt T f t 1 1 1 t 1 1 1 1 e 
ES TKK) COG KKK, Lec CK EK EKOR. SK! KUM o 
f T 1 1 1 

NO KOXK Kl, KIK, +K] KI, KE 


1 1 t t ! ! t 
KU K+K.K, I+ (KI+K3) Kyt (K1tK3) Ky e e (KIEK) Kat (K+K) KL 





1 + 1 1 f t 1 
KIK +K K, KPKK +(K +K DKI n. (K+K +K) K + 
: (K+K +K) Ki 
1 1 t 1 1 1 1 t 
ER ARO HK KY, (RRA 
T 
EA 1 A 


4 2 





ic, KK 


152 Sify o EK 





lnn 
f En 
KR, (K+K) Ko E (Ki KK. 
= AS E +K + 
KK, (K\+K,IK, (K K, K.IK, ze K3 KK 
Ya T att ge eis 
KK, (K+K) Ky (K+K,+K.IK, ; 2e Ka ei 
and 
eT 
E, 
= + 1 1 t 
Y - LK,» Kyo Ko K J E, LK, Ky > Ke 
En/2 
where 
5 Bentl * (n/2-1)m*1 
ic E. *m*2 ES een? 
= = E E 
El Eo 5 
ml “om _| =a im 
2 > > > x r = 
T Naentity matyıx, and Ey > E, Bayo? Er > E,» 
have the same dimensions as E E | 


PAZ 


and m corresponds to number of inputs or outputs. 
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ue ntescestinpso Note that the simplification of the 
state equations can be carried out by partitioning the matrix, 
or by only keeping a part of the original matrix as a simpli- 
fied model. In other words, if a simplified model of a two 
dimensional matrix is required which is equivalent to a 
2 x (2m) order transfer function, the upper left nand corner 


of the original matrix is taken as the simplified model. 





Therefore: 
1 t t t 1 1 1 + t ^ 
KK K3 Ky ke 1+X4K,+K,K, K Ky K- Ky Er 
ty! t ty t t t t t = 
KK, (K +K) K RU KotK Ky I+ (K+K JK t (K+K) K, E, 
IK KK E I 
Nm lu en U 
DEUS CK. *K42K E D 
i 2 1 3° Wu 2 (80) 
Ei E 
= t t 
Y = LK», K, J TON. Kid 
Bo Bo 
Coa) 
where the vector Er ED Er > E, Er E, keeping the same 


dimensions as were given for the equations (78) and (79). 
The method used for the simplifications of a transfer func- 
tion and state space formulation based on the Cauer mixed 
form holds in the same fashion for the Cauer first and second 
forms. 

For effects of simulation Figure 7 shows a diagram for 
two inputs two outputs fourth order system for the Cauer 


l3 xed Form, ın samılar Form for the Cauer first and 


"m 
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PAZ 
Figure 7. Simulation diagram for fourth order system with m = 2 
for Cauer matrix mixed form. 
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second form a simulation diagram can be obtained as given 
inerneures 8 and 9. 

hes important to notiee the number of reductions, ca, 
which one can realize are constrained by the number of inputs, 
m, and the order of the system, n, where q = n-km for 
E 2... and q > 0, 
ERSENUMERICAL EXAMPLES AND GRAPHICAL COMPARISON OF THE 

BHRES CAUER FORMS 

iin example 2.1 

Consider two inputs, two outputs fourth order system 

where a reduced second order system is required. 


The state space equations for the system are: 


0 0 95.2 5 = 1121225 lo 0 
0 0 2225 =), 125 0 0 
x = wur u 
-1 DES -1.5 -0.5 Der 
=].5 =0.5 -1.0 -4.0 5 05 5 
l P (82) 
i 2i Au l 
Y = x 
0 jt -1 d 
(83) 


The transfer function for this system is given by equation (84). 
-l 
stl s-1 25* Ta Sur? s-1 


iM sz 2 
-2 stl -4s - 4s - 1 2s +65 +3 


(84) 
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Figure 8. Simulation diagram for fourth order system with m = 2 
for Cauer matrix first form. 
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»imulation diagram for fourth order system with m 
for Cauer matrix second form. 
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From the transfer function, the matrix Routh's Array and 
Hae following matrix quotients for the three Cauer forms are 
obtained. 

Therefore, reduced second order models are: 


Gaver first form 


ME 5 -il , 
eS y! ES 
1 m 3 2 1 5| 13 


Cauer second form 
DES 0 
O ARAS 5 


2 i 
H, = H 
1 1 4 2 
1 = 
-1 S 


Cauer mixed form 
where their respective transfer functions are given in equa- 


A 
Ha 
pl 
| i | 
= N 
pt 
A 
pa 
n 


pens (84), (85), and (86). 


Sauer first form: 


Sce b coc 05 1 
Fa] (S? - w 
SROS s + 2.5 zo pgs is 


(84) 


Cauer second form: 


us + 2 -l 1 
Poo?) = EE e 
xs I ASA ds F 3s T 5 
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(85) 





Cauer mixed form 


IS? Ss = l 1 
ER (5) = > 
s + 1 Se re UE. Shot 


Figure 10 presents the step responses of reduced models just 


(86) 


developed, where (10a) corresponds to the output (1) and (10b) 
eerresponds to the output (2). 
ae exemple 2.2 

The following example shows the reduced systems in 
Cauer first and second forms are only stable for a second 
order approximation, their fourth order approximation pro- 
vides an unstable response. Whereas, the fourth order 
approximation for the Cauer mixed form is completely stable. 
Furthermore, the Cauer first form gives a poor approximation 
in the steady state portion as is expected. Consider the 
following sixth order system with two inputs two outputs as 


given in equations (87) and (88). 


0 0 20 0 0 0 0 0 
0 0 0 M) 0 0 0 0 
x =| 0 0 0 0 des) 0 > O DE us 
0 0 Š 0 0 0 Leo 0 0 
zip 9.96 -108.5 19.6 -]6.5 Hu 2 -l 
49 -19.96 107 -39.6 18 -10 -2 2 
l (87) 
HO Dal RE 0) 0:52 
Y = > 
1.0 4 2 I p 
(88) 


50 








* SWAOJ Jane) 
39442 34] UL U9ISÁS ABPAO puooss paonpau pue WOISÄS Uapuo YIANOY [PULÉLAO. 3YY 404 S3nd3nQ "OL SANBLY 











(9) Ub) 
IVES Leet en 0 71830) 1 ON 1094 SN IEC la rG o -i^ 
RIVAS Ia Y 2 ez SR i “ON NY ADA A 3C > 3% ` 
ei. Tt » 5 ^0 à co” fac? 00 4 ooe C Oui Oeo e a a eee 
: | 
m 
2 = 
N 
| 
o i a 
= J> 
E . 
t 3 - 
j e A LL — 
WO SALJ Aone 5 
ne J a, WJWOJ puodas Jane) —» i 














1W9ISÁS eurbrao ——> : K 6 
10) paxjw 4one) —Z oe oe 
Be EE Ei e | 
UMOJ puooes Jane? : 
: WAOJ 1SJ1j Jane) ] | 
a | : 
: n 
o i 3 
E 
Lə 


EN 





Miertranster funetion of this system is given by equation (89). 


S t+ 10 0.282 AS O 


Ts) = ) 5 
SE Zas T SUL 


x t De AO O SS Be + Bao: TZ -1 


A a a 0 


(89) 
or if the inverse is to performed 
5 y 3 2 
-0.25s + 0.ls + 14.95” + 79.25s + 209.755 + 99.95 
mcs) = 
> T es + CE t NE EG 4039.55. 778 
DES E SE + dom E NN DX CES p H 
ENS - E M b 93.588“ + 2uu.76s + 99.98 
en 
6 5 T 3 | 2 
EET: 173-25s + 130.558 + 562.535 + 1349.13s +1049.875+2,9,99 


(90) 
a. Cauer Matrix First Form 


Performing the respective Routh's Array: 


1 0.5 10 1 55 0.2 

1 1 3.5 6 ZO 
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um tf = 
EL Ñ n 
5 1 0.197 
25556 SS 02.076 
m m 
Hh = A, 
1996 3326 -0.28 
1.64 =D a, 0.29 
TEE lo «« 
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The state space representation in 


given by equations (91) and (92). 
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respective H'i's are: 
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The reduced fourth and second order model are obtained by 
partition of the respective matrices in equations (91) and 
(92). For this case, the fourth order model is unstable and 
the second order model gives an untolerable error in the 
steady state portion of the response, so this method is not 
applicable for this example. 

Figure 11 shows the step response for the original 
system and the reduced second order model. 

bee Caller Matrix Second form 


Performing the respective Routh's Array: 
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and from the Routh's 
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The state space representation in 


given by equations (93) and (94). 
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Cauer second form 
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The reduced fourth and second order model are 


OPOPOP 
HOHORO 


ss) 


(94) 


obtained by partition of the respective matrices in equations 


es) and (94). 


moumehn order model 15 unstable and the second order model 


offers a good approximation to the step response of the 


S 


For this case like the Cauer first form, the 








original system. Figure 12 shows the outputs of the orıginal 
and reduced second order model for a step input. 
c. Cauer Matrix Mixed Form 


Performing the respective Routh's Array: 


D» 55 (07 
BS LO S 
2L i Or. 
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-54.31 Es S 5 4.04 
21-837 20.9 
-7.2 -0.81 


and from the Routh's array the respective Ki's and K'i's are: 
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hod fog fs. 
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O O 
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2.505 ERG 2.5 1 
ER t = 
py = Ss 
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€ t u 
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2.24 1.079 TEITCEBENGI 
10.738 -6.428 DESC qe SIS 
= 1 = 
K3 = K3 ~ 
-27.468 14.532 -6.77 2.556 


The transfer function for the reduced fourth order system 
is given in equation (95). For this model the step response 
is completely stable and the approximation has an error of 


only 0.5% over the entire response of the system. Figure 13 
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shows the outputs of the original and reduced fourth order 


model for a step input. 


-0.012s° - 0.063s* DES re a 


fics) = 3 2 
ebs rer 7325 vt O0 273s F 0O03 
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O SAR 


a 


0.7955 + 8.215s% + 8,7188 + 2.431 


60 


(959 





ime DERIVATION OF LOWER ORDER SYSTEN FOR THE LINEAR 


i SA EE ee EE 


REGULATOR PROBLEM VIA CAUER FORM 


Pee SbINEAR REGULATOR SYSTEM 
R. E. Kalman [12] has shown that when a system is 
described by equations ( 96) with a performance function 


given by equation (2397) 
i= Ax + Bu ( 96 ) 


5 [xl DESEE UE R u] dt (097. 


C4 
N 


and if the plant is completely controllable then an optimal 


u* exists and is given by equation (98 ) 


F 


n oz Kx ( 98 ) 


where the feedback matrix X for the controller depicted in 
Figure l^ is a constant matrix as t-. 

In recent papers, Goldman [10] and Aoki [18] have 
derived a way to obtain near optimal solutions using the 
"aggregationmatrix", to transfer an arbitrary state x to the 
origin of the state space while minimizing the criterion 
ae tion given in equation ( 97 ). 

The present work will show that the reduced order optimal 
regulator can be obtained where the original system Is 
translated in Cauer form and also a near optimal solution 


can be found through a mere partition of its optimal solution. 


6l 








| REGULATOR, ___| 


Figure 14. Linear Regulator System 


feeds STATE VARIABLE REPRESENTATION BY CAUER 
SECOND FORM AND ITS OPTIMAL SOLUTION 


It is known that the standard state variable representa- 
tion of a linear time-invariant system is given by equations 


Woo) and (100). 
x = Ax + Bu (33) 


|j wer (e 


where x is a column vector nxl given by LX, > fund E A 


is an nxn constant coefficients matrix, x 1s a column vector 
T : ES 
Mx] given by Lx,» X55. X] , B is an nxl constant coefficients 


02 





Metrix, u is the input, y is the output and C is 1xn constant 
coefficients matrix. 

Goldman (10), Chin and Shieh (11) have proven that the 
system given by equations (99 ) and (100) can be represented 


in Cauer second form as shown in equation (101) and (102). 


z= Fz + Gu (191%) 

y = Mz (102) 
where, 

Fs pap. (103) 

EE PB (104 ) 

M= cp”! (105) 


where the matrix P is an nxn upper triangular matrix and the 
elements in the triangle are copied directly from the elements 
of the Routh's array, where the nth row of the P matrix 

is the (2n+1)th row of the Routh's array and the states 


variables x and z are related by equation (106). 

Z = PX (106) 
lune system in (99 ) and (100) is put in transfer function 
notation the result is given by (197). 


7(s)@21ta2284a238 ra2US 4. rang" 


> 3 = NOS 
alltal2S+al3S+taluS”+...+a2n+1lS 
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The Routh's array can be formed as 


eed? ars ra. ‘alm "alintl 


eee a2 a3 Ce. Salen (108) 


and the elements of third, fourth, and subsequent rows can 


be evaluated from the following algorithm. 


ajk = aj-2, k*1-fj-2aj-1, k-*1 Wace ee eh 
E Let: 
(109) 
E = 
fp = === Delos m 
SUN (11:07) 


The elements of the matrix F can be obtained also from 
equation (111) where ap,l and ap*l,l correspond to the 
elements of the first column of the Routh's array and then 


F is formed as shown in equation (111). 


wor] PTT Fo AH TEE 
ITI TA(ELFES) fSCTITÍ3) > seems LTS) 
EA COII VCEHCEITEISO S fOCEITIÍSGTEDO...fnCflitt3trt5) 


SU or SR ES nd LAS hs n=) EN 
( ) 


Ihe elements of the matrix G correspond to a column vector 
of nxl of 1's and the elements of the matrix M can be formed 


meso by ip's where p=2,4,6,...2n as shown in equation (112). 


M = [f2, fh, f6,...f2n] 
(112) 
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The performance function given in equation (97 ) can be 
translated to Cauer form by using the relationship given 


iemeguation (106), i.e., 


nl 


SCs Q zt + uti RU+ Jat o 


Ci 
H 


(114) 


oO 
tt 
"U 
oO 
"I 


u® = TS Ze CON] a) 
where 
-] 
E Ton 
Ko K P (116) 


So the equivalent system in Cauer form given in Figure 1 


can be represented by Figure 15. 





REGULATOR IN CAUER _ l 


Figure 15. Linear Regulator in Cauer Form. 
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Note that K_ is the feedback gain matrix of the Cauer 
system and u* is the optimal control law of both systems 
(Figure 14 and 15). 

From the theory of optimal control, z* and u* are found 


by solving the set of necessary conditions: 


Zu = E MU put) 
dp* 
ee (z*, u*, p*, t) 
gz 
Torvall 
iE EE 
0 are (A ER pr, t) 
2 5 
nz = Gal. už, p*, t) 
(117) 
Where the function H is called the Hamiltonion and is 
defined by: 
> A ~ 
H(2,u,p,t)=slz Q z*u RU]*p. [Fz*Gu] 
QUES) 


and D 1s the Lagrange multiplier or Costate state. 
Substituting equation (118) into equation (117), the 


canonical system in equation (119) is formed. 


E „IF exe || as 
$^ -Q, -F p^ (119) 
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Fortunately, for the Optimal Regulator problem, it is not 
necessary to solve these equations. Kalman [12], Tyler, 
u.S. and Tuteur, F.B. [19] have shown that when the optimal 
@entrol u™ is generated by equation (115), the solution 

LOP Ko is obtained from equation (120) where the matrix, T, 


is the solution to equation (121) in steady state. 


MWER CT (120) 


i 


F°T+TF-TGR + 


© 
u 


G^ TQ, (121) 


Equation (121) is the steady state form of the Ricatti 
equation. 

For the single input case the weighting matrix, R, is 
a scalar. When tnis is the case, a lower order optimal 
linear regulator can be found just by mere partition of 
Eo. M, Qc and K matrices as shown in equation (122). 


Hence the new system is depicted in Figure 16. 





Figure 16. Reduced Linear Regulator in Cauer Form 
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E 
M= [f2 £u ı £6...f2n] IESU RD eun 
i 
Ver Send 609 4o2n 
K_= [k K A IP c 
e eu C12 ein |Sen1 don? N 


(422) 


The reduction scheme outline above has left u and uns 
unchanged. Therefore, gus is tee suboptimal Solution. For 
all three systems, that is the original system, the trans- 


formed system and the reduced system. 


EE REDUCED ORDER SOLUTION 


Based on equation (122), the reduced system is given by 


zr = frz + Cru (07:35 


y Mrz 0424) 


where Fr, Gr and Mr are the partitions of the matrices F, G 


and M respectively. 
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Goldman [10] has shown that the original system is 


related to the reduced order model by the following equations: 


a + 
Were the matrices € and € 


Ar 


Ar 


+ 
eÄe zr + eBu 


er 


EX 


Zr 


Ar 


ame poa iton ot the Pande? 


GEEZ So 


(122065) 


(12:70) 


(128 ) 


matrices in rectangular form respectively, such that e is of 


to. 
order (rxn) and e is (nxr), where, n, is the order of the 


original system and, r, is the order of the reduced system. 


then 


oF. 


Note that from equations (98 ) and (06) 


Or 


From equation (127) and ( 98) it is easy to show that 


Kr 


= 
= 


it 


Ke” 


Ket 


Kre 


KP 


aul 


Zl 
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(129) 
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Then from equations (130) and (133) 


Nm K Ir (134) 
where 
+ 
Ber =- Pes (135) 


The matrix Ir has dimension (nxr) and is a special type of 
identity matrix where the last (n-r) rows are zeros and the 
ASES xr components correspond to the identity matrix, In 
the same way, it is possible to show that the different 
Relationships that exist for the weighting matrices, Q, Qc 


Feen are given in equations (136), (137), (138) and (139). 


pop Gee) 
Qr = e Qe Grea 
a e P OcPe (138) 
Qr = Ir? Dc D 


1. Example 3.1 
The transfer function of a system is given in equation 
(140). It is desired to control the system in such a way that 
a performance function given in equation (141) is minimized. 
Due to the systems complexity, a near optimal solution is 


desired. 


H(s) = mn (140) 
Sales 
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T = Peix-OxtU RUIdt (141) 


where 
Del 
Q= JO 4 OQ Rea? 
DTO l 


0 1 0 0 
A= JO 0 aL EUN) eco E | 
UI 5. «5.1 dy 


The Routh's array 


SP Br ALSO 
D 0 Du DO 
p qu 
-1.09 -0.179 
PEN 
9.932 
HO, 
From the Routh's array 
55b c 5 au 
P = 0 AZ 120 
0 0 jd (142) 


a 





20-2107 069032 
02195 0.192 
Teu (143) 


then 


a e 
2 
€ S 
Ss? 1080 (144) 
.179 


=0.2102 


0193 


M 
i 


(145) 


The Qc matrix obtained from equation (136) is 


0.159 -0.188 ! 0.029 
l 
eS -0.188 0.372 | -0.183 
Ecce LULA a 
0.029 -0.183 1.154 (196) 


The Fe, Ge and Me matrices from equations (103), (104) and 





(105) are 
0.105 ! 0.015 
[ 
Foz ESSE S 
MEAE Mem aeree im eod 
Bass Ed GG 
D 
Ges | 1 Me 178620. 21120.032] 
--- 


1? 





The reduced system is given by equation (148) and the 


corresponding matrices that describe the system in equation 


ES). 
CONS LO 
s*+1.065s+0.0965 (T33) 
S 0.105 
FP: 
-0.089 reo 78 
17 


Gr= | Mr- [0.1786-0.211] 
a 


DIES -0.188 
-0.188 0e 72 (349) 


The computer solution for the three systems for the feedback 


gain are shown in Table I. 


ABEL 












OPTIMAL VrEEDEBACK MATRIX 


227337293425 2.08.7256 0.41034 


ESUER FORM reo 3919509 D-20250 S33 H0. 065519992 
REDUCED CAUER FORM Dogs ZH 0.02288695| | 


Sr SE 
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Table II shows values found from the reduced system 
back to Cauer form and original system with the equation 


a and (132). 


TAB ees 


SYSTEMS NEAR OPTIMAL FEEDBACK MATRIX 
12273008 PS BS 


Dco db S 0:02 2506 


REDUCED CAUER FORM 0231965 02.022580 


Note that a near optimal feedback matrix is obtained directly 





from equation (131) but not to the Cauer system since the 
reduction e discarded the last terms. 

Figure 17 shows the step response for the optimal and 
reduced linear regulator. Figure 13 shows the optimal and 
suboptimal control laws. Figure 19 shows the step response 
for the optimal and suboptimal linear regulator. 

mee Example 3.2 

The methodology just developed is applied to a simple 
model of a nuclear reactor power generator [20, 21]. 

The heat generating process of a nuclear reactor is 
dependent upon the mechanism called fission (a fragmentation 
of matter). The power generated by this process is directly 
related to the population of neutrons, n(t) and can be 
described by the following differential equation (developed 


from a diffusion balance equation). 
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5. 17 uode SA o UNS E 

meiMAe AND REDUCED LINEAR REGULATOR 

d 

S|/ 

00 2.00 4.00 6.00 8.00 102 00 12.00 14.00 16.00 15.00 
XSCALE= 2.00 UNITS/INCH RUN NO. | 
“Seelen af UNITS/INCH PLOT NO. 1 
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SEE EL EN 
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dg NDS 


18.00 





305 ea) aud etus (150) 

Eee) e (nt) NEC Gis) 
where 

pee) = 6 Ke(t)=an(t) US 


The variable 6k(t) is the input to the process and is given 
“me name "reactivity". It is clear by inspection that 
§k(t)<B for stable system (in a linear sense). The 
variable c(t) is a measure of the concentration of fragments 
(precursors) that produce delayed neutrons according to a 
time delay (1/A) called the "half-life" of the precursor. 
iimemamput, 6ke(t), is the control (reactivity) that is 

Ir @ejlated with the control rod position and, d, is a 
temperature feedback (reactivity) coefficient. 


The parameters for this system are: 


"NEEDED 
= ine sec. 
B = 0.0065 


A = 0.1 sec. Tt 
at t=0 n(0)=10kw (as the operating output in steady 
state Conditions) 
The control problem is stated as: 
Find the Optimal Control Policy u*(t) that will 


transfer the power level n(t) from the operating level 
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n(0) = 10kw to a new level n(t) = 50 kw where 
u(t) 2 8 ke(t) Gaos 


using the performance measure given in equation ( 97) where 


des 0 0 
2-10 0 0 
0 0 25 Gist 


and compare this solution to the near optimal one obtained 
from the reduced order model. 


From the state of the problem the following diagram 










Poe drawn 
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i KINETICS 
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=K E 
L 
Figure 20. Block Diagram Optimal Control of a Nuclear 


Reactor. 
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The double subscript, ss, means steady state and the 


new definition of state variables are: 


pa) 
DT) 
po (tc) 
The easy. 
n(t) = 
c(t) = 


Sk(t) = 


Ae initial 


from the state variables definitions, equations (155), 


(156) and (157) and equations (150), 


Rt) 


c(0) 


EOD) = 


me (0) = 


ez) a= 


Sk(t) = 


6k(0) = 


state conditions are 


0 


ne) 
aCe) 


6ke(t) - Ske. 


0 
0 


condi ions at Nn t) 


È 
1 


> 


1 


640 


c(0) 


-2960 


an(t) 


anlo) 


C and Lolo 
W eda a) 


B m 0) 


AN co 


(T595) 


(156) 


Cio) 


(1589 


525 
(160) 


CLOT) 


(le2) 


(163) 


(164) 


(165) 


(166) 


(16%) 





ak(0) = 10 
“ss 
QUSS = dno. - Àl n r B 
SS 
8k E 
SS 
BAD) s ó5k(0) -Gk 
Ss 
_ -4 
Zo (0) = -4x10 
x1(0) = n(0)-n 
ss 
mi O) = =40 
then 
[x1(0) -40 
a0) = x2(0) = -2560 
x3(0) ln 


and the optimal control law U*(t) will be 


Mee Ct) = -Kx(t) 


same as given in equation (98 ). The system in variable 


form is given by equations (99 ) and (100). 


(160) 


(205199 


(MAGO) 


(1/2) 


CI 722) 


GE 


(174) 


WES 


Since equations 


(150) through (157) are non-linear linearization is required 


(see appendix A). The linearized system is represented by 





1 SS 
m CAU 0 
(en ¿+8) X 1 
s E 8/1 -À 0 x + 0 u 
0 0 0 ak 


: 
| 


ye) all 0 dC E) 


where 
y 
-6.9 ue 5x10 
ABB e 0 
0 0 0 
BS 


Solving the Ricatti equations, the matrix K is given by 


Ccemarion: (181). 


L 6 


Kae eS ie 8x10 7.05] 


The transfer function for the optimal system will be: 


Blin Corot 


3 2 


H(s) = — > — 
b Pistol costo 
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( 178) 


(179) 


(1909 
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The Routh's array for this system is: 


ny AWS 61.9 Ors 1.0 
5000 50000 0.0 
0.3 14.05 120) 
46539 -246.3 
n r5 720 
-3533.58 
Eo) 


Prom the Routh's array 


20.3 1005 1.0 
P= 19.0 ido ae 1.0 
DONEC SOS 1.0 (183) 
0.049 -0.049 -0.0004 
pt = o.o 0.07 ORO 
050 Du) O (184 ) 


and the rectangular partition of P and ae matrices give 
2023 14.05 Tz 0 
€ = 
0.0 IS 1.0 (135 ) 
0.049 -0.049 
ce = [0.0 0.07 
0.0 0.0 (186 ) 
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and from equations (125) and (126) the reduced second order 


suboptimal system will be 


ä -0.205 =2.135 a 
zr(t) Zr) mE 


-0.205 -4.169 l| 
(187) 
wee) = [246.305 SSR ZnT) (188 ) 
The corresponding transfer function is given in equation 
(189). 
SS 
CUI (189 ) 


In Figure 21 is shown the step response for the original 


optimal and reduced suboptimal system. 
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HS. 21 EXRMPLE 3.7 - STE 
Bem OINAL AND REDUCED OPTIM 
S 
S 
S 
S 
S 
m f.original system 
reduced system 
a 
S 
S 
S 
c5. 00 8.00 16.00 24.00 32.00 40.00 
XSCALE= 8.00 UNITS/INCH 
YSCALE= 4.00 UNITS/INCH 


09 


48.00 


56.00 6 


RUN NO. 
PLOT NO. 


. 00 


1 
1 


12.00 





IV. CONCLUSIONS AND RECOMMENDATIONS 


Comparative analysis of Cauer forms methodology for 
multivariable system reduction is established. The developed 
methodology is based on the Cauer matrix generalized form, 
which offers the closest approximation to the original 
system. The proposed methodology and because of the nature 
of the Cauer mixed form shows it to be superior to any other 
method proposed to date since it provides satisfactory 
results for both the transient and the steady state portion 
of the system response. The methodologies in state space 
as well as in the S domain are developed. For the basis 
of comparison, reduced order models using the three Cauers 
forms are obtained for two different examples. The results 
clearly show the superiority of the Cauer mixed form over 
the entire frequency range of system's response. The 
proposed methodology is algorithmic therefore, it is 
amenable to digital computation. 

A lower order optimal linear regulator can be obtained 
by mere translation of the original system to Cauer second 
form and their partition of the different matrices as shown 
in Section III-B. A suboptimal feedback matrix for the 
original system can be obtained by multiplication of the 


lower feedback optimal matrix by the rectangular partition 
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Stecne Matrix P. The responses of both systems, the original 
and the reduced, as well as optimal and sub-optimal systems 
are extremely close. 

The results presented here are encouraging, there is 
a distinct need for future research, in particular in the 
reduction of systems with any given number of inputs or 


SREE DUT S . 


o 





APPENDIX A 


Linearization of Multivariable Systems 


Let equations (A-1) and (A-2) be the state variable 


modeling 
pene (ee Uh) (A-1) 


Eau) (A-2) 


< 
H 


For the system described by equations (A-1) and (A-2) is 
operating at steady state conditions(constant input Use» 
producing constant state ie and constant output Ve E The 


combination of these produces 


DE c Uu i2 (A-3) 
ss ss 


le) (A-4) 


g(x E 


Yes SS 


If the system is perturbed by either drawing the states or 


the inputs, the system motion satisfies 


DE ce toe O t5 u) (A-5) 
Ss ss 


y _Hy 


ac tur Bm) (A-6) 
ss ss Ss 


Both functions, F and g can be expanded ina Taylor series 


about the points x ao aed resulting in the following 


representation of the system equations. 
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e i OD E U s? tA x186 uta (6 x ,ó u) 


(A-7) 
UE o y-gix... u$, *C8 x*D6 u*8 (6 x ,$ u) 
(A-8) 
where: 
= nxn matrix 
B = nxm matrix 
"Ec vx matrix 
Income matrex 


n = order of the system, m = number on inputs, v=number of 
mucputs. 

The functions a(óáx,du) and 8 (6 x,6u) represents all 
second order and higher terms in the Taylor series expansion. 
Substitution of equations (A-1) and (A-2) evaluated at 
vu BO and neglecting second order terms and higher 


S9 S 


terms, yields the following perturbed equations of motion 
6x - Ax + Bu (A-9) 


cx + Du (A-10) 


Ó y 


The two equations above approximate the dynamic behavior 


of the system about the operating point en Ute The 


elements of the motion A, B, C and D are given by 


Se o 1.2... Ede o cala 


eg 








IR 


Bij 


I 


er] 


il? 
© 

a 
H 


Dij 


Ihe representation of equations (A-9) and (A-10) is given 


by the following block diagram. 





Figure A-1. Block Diagram of a Linear Model. 
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